Wave-packet trains of a time-dependent harmonic oscillator 
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By using a test-function method, we construct n exact solutions of a quantum harmonic oscillator 
with a time-dependent "spring constant" . Any n-th solution describes a wave-packet train consisting 
of n + 1 packets. Its center oscillates like the classically harmonic oscillator with variable frequency, 
and width and highness of each packet change simultaneously. When the deformation is small, it 
behaves like a soliton train, and the large deformation is identified with collapse and revival of the 
wave-packet train. 
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Classical and quantum motions of the time-dependent 
harmonic oscillator are practically important, because of 
their physical realization in a Paul trap Q, 0. For 
some periodic and time-dependent "spring constants", 
the classical motions are governed by the well-known 
Mathieu equation 0. Most of works for the system 
are of classical dynamics, such as the stability analysis 
P|, order and chaos 0)0) an d so on - I n the ultracold 
ion case, quantum motion becomes important for prac- 
tical applications 0, [Hj- For convenience' sake, many 
works based on quantum mechanics, say, the quantum 
computations [ll|, and state preparations 0], E3) 
were performed by adopting the pseudopotential model 
@ with constant trapped frequencies. There exist a 
few of works on the quantum motions in a Paul trap 
with periodically variable trapped frequencies. Cook and 
coworkers proximately treated the time-dependent quan- 
tum oscillator by perturbing about a harmonic oscillator 
solution for a time-independent, effective potential [l5| . 
Combescure, Brown and Feng gave similar ex- 

act solutions of the time-dependent system by employing 
different methods respectively. Their results showed that 
quantum-mechanical solutions of the system depend on 
the corresponding classical ones. 

Very recently, Moya-Cessa and Guasti reported a co- 
herent state solution of the time-dependent quantum har- 
monic oscillator pjj that may be associated with sum- 
ming Brown's n exact solutions. For a time-independent 
trapped frequency, we suggest a test-function method to 
get n exact solutions of the system |2fJ. In the present 
paper, we shall adopt this method to find n exact solu- 
tions of the time-dependent quantum harmonic oscilla- 
tor. These solutions describe n wave-packet trains that 
propagate and breathe simultaneously. By selecting ini- 
tial conditions to fix centers of the trains, the result is 



reduced to Brown's one [lTj. And the ground state with 
n = may be related to Cessa's coherent state |l9). 

We consider a single ion confined in a Paul trap that 
provides an oscillating quadrupole potential with the re- 
sulting quantum motion governed by the time-dependent 
Schrodinger equation 
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where the "spring constant" is a periodic function of 
time, k(t) = U 2 + V cos uit with ui being the rf-driven fre- 
quency and U 2 , V the trap parameters. We have adopted 
the unit with m = h = 1 and will normalize time in unit 
2/u) so that k(t) = U 2 + I^cos2i. The spatial coordinate 
x and probability density |^»| are normalized by the har- 
monic oscillator length 1^ = \J H/Jmuj) and inverse one 



of Eq. (1) in the form 



h , respectively. Consider the test function as a solution 



ipn = a n (t)H n (£) cxp[b(t)x - c(t)x 2 
£ = e(t)x - f(t). 



f(t)/2] 



(2) 



Here a n (t), b(t), c(t) are the complex functions of time 
and e(t), f(t) the real functions. Applying Eq. (2) to 
Eq. (1), we arrive at the equation 
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H n = 0. 



(3) 



Noticing the Hermitian equation d 2 H n /d^ 2 
2£dH n /d£ + 2nH n = 0, Eq. (3) implies 
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ic = 2c 2 - k(t)/2, ib = 2bc, 
ie = 2ce — e 3 , if = be — e 2 /, 
ia n /a n = iff — b 2 /2 + c + ne 2 . 



(4) 
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The first of Eq. (4) is a complex Riccati equation, 
which can be changed into a complex equation of a clas- 
sical harmonic oscillator 

Cp = -k(t)<p= -(U 2 + Vcos2t)ip, (5) 

through the function transformation c = tp/(2iip). This 
equation is just the Mathieu's one whose exact solution 
can be taken as the form of infinite trigonometrical series 



Q . Generally, for given physical parameters U 2 , V, u> and 
initial conditions to determine coefficients of the infinite 
series is difficult. In order to overcome this difficulty, we 
here shall employ an iteration method to seek the series 
solution, through an integration equation. Setting real 
part and imaginary part of ip as ipi and <p2 respectively, 
they satisfy same equation (5). The corresponding inte- 
gration equations read as [21| 



<pi = Acos(Ut + a) + VU 1 ^sint/t J cos Ut cos 2ttpidt — cos t/i J sinUtcos2t(pidtj , 

ip 2 = Bcos(Ut + (3) + VU~ 1 ( y sinUt J cosllt cos 2tip 2 dt - cos Ut J sint/t cos2t(p 2 dt] 

I 



(6) 



where A, B, a and /? are some real constants associ- 
ated with initial conditions of the classical oscillator. In 
quantum- mechanical treatment, there exist many possi- 
bilities on these constants that correspond to different 
states ip. The two equations in Eq. (6) can be directly 
proved by inserting them into Eq. (5). Applying Eq. (6), 
we construct the solution of Eq. (5) as 



P(t) 



<Pl + W2 = P(t) 



9(t) = arctan 



f2_ 



(7) 



The real functions pit) and 9{t) will be use as the known 
functions in follows. So far, most of experimental works 
and theoretical analyses on the Paul-trapped ions are fo- 
cused on the first stability region where the trap 
parameters U 2 and V are small, U 2 < 1, V < 1 and 
V ~ U 2 < U . In such case, we can treat the terms pro- 
portional to V of Eq. (6) as perturbations and adopt the 
iteration method to produce the perturbed solution. 

Returning to the transformation between ip and c 
yields 



2itp 



■ P_ 
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(8) 



Substitution of Eq. (7) into Eq. (5) yields equations of 
the phase 9 and module p as 



-26p/p, p = P 2 - k(t)p 



with the first integration 

c = P 2 e = 



(9) 



(10) 



Combining Eqs. (8) with Eq. (4) and applying the rela- 
tion (10), we easily obtain 
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VP I L\ 2 / 4c 
where bo and Aq are arbitrary constants. Insert- 
ing these into Eq. (2) leads to the exact solution 
tp n (x,t), and the normalization condition J \ip n \ dx — 

A 2 ^ jH 2 (£)eM~e)dt - A 2 ^c^2 n n\ = 1 gives 

the constant Aq = L/co/ '(v / 7r2"n!)] 1 / 2 such that from 
Eqs. (2), (8) and (11) we have the normalized wave- 
function 



4>n 
Rn 

8„ 



R n (x,t) exp[iQ n (x,t)}, n = 0, 1, 2, ■ 

1/2 



V¥2"n!p(t) . 
p(t)x 2 box 
2p(t) ~~p~(t) 



Bin ^t) + sin[20(t)l -(1+ n)9(t). 
4co V 2 / 



cos 9{t), 



(12) 



The module R n of exact solution (12) for different quantum number n describes the wave-packet trains con- 
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sisting of n+1 packets. Orbit of center of the wave-packet 
trains x c (t) is given by £ = 0. Applying Eq. (7), from 
£ = we have the orbit 



— p(t) cos 9(t) = —(pi, 

Co Cq 



(13) 



which is proportional to real part of the complex solution 
(7) , namely the orbit of a time-dependent classical oscil- 
lator. Comparing our Eq. (12) with Eq. (24) of Brown's 
paper [17]. we find that if one let the arbitrary constant 
bo be zero, our solution agrees with Brown one. The 
orbit equation (13) indicates difference between behav- 
iors of the two solution: our wave-packet trains oscillate 
their centers, but the centers of Brown's wave-packets arc 
rested. 

The function p(t) appearing in £ describes widths of 
the wave-packet trains and each packet. For any train 
the average width of the packets is p{t) / y/co. Same func- 
tion appearing in radical of Eq. (12) governs highnesses 
of every packets. So we call p(t) the function of width and 
highness. When the changes of the widths and highnesses 
are small, behavior of the wave-packets seems to be that 
of the soliton trains. The greatly variable widths and 
highnesses show collapse and revival of the wave-packet 
trains, like behaviors of multiple breathers. The normal- 
ization condition implies that the broader wave-packet 
train is associated with smaller mean highness and the 
narrower wave-packet train corresponds to larger mean 
highness. 

Phase of the exact solution (12) contains the term 
(1/2 + n)8(t) that infers average energy E n (t) = 
(ip n \id / dt\ip) to be proportional to (1/2 + n). In the 
processes of motions, the wave-packet trains may spon- 
taneously transit from states of higher average energies to 
that of lower ones. Some perturbations also could cause 
transitions between the states of different quantum num- 
bers. The ground state ij>n of Eq. (12) may be similar 
to Cessa's coherent state [jjj. We shall numerically il- 
lustrate the analytical results for the system parameters 
U = 0.5, V = 0.05 as follows. 

The solitonlike trains To show the solitonlike be- 
havior, we require larger amplitude of Eq. (13) and 
keep small change of module p such that the wave-packet 
trains described by Eq. (12) can propagate through 
greater distance with small deformation. For simplicity, 
we take the parameter set a = 0, j3 = —tt/2, A = B = 
Co = 1, 6o = — 10 and make iteration from Eq. (6) only 

to the first order. That is, let tp^ and < / 9 2 °' > be cos Ut and 
sin Ut respectively, and use them instead of tp\ and (f2 m 
the corresponding integrations of Eq. (6), obtaining the 
approximate solutions 

(fi w cos 0.5t + 0.05 sin 0.5i sin 2t 

+ 0.025 cos0.5i (1 -cos2t), 
tp 2 w sin0.5t + 0.025 sin 0.5i (1 -cos2t). (14) 

Applying the first of Eq. (14) to Eq. (13) yields the orbit 
of wave-packets' center. Obviously, the center of wave- 
packet train oscillates in the spatial region — 10 < x c < 10 



with period Att. Combining Eq. (14) with Eq. (7), we 
obtain the function of width and highness pit) whose 
time evolution is plotted as Fig. 1. In the considered 
case the module oscillates with very small amplitude that 
implies width and highness of the wave-packet train being 
changed also quite small. 
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FIG. 1: Time evolution of the module pit) of the classical 
solution that exhibits small change of the solitonic width and 
highness with period 2n. The space-time variables x c and t 
are normalized in units lh and lj -1 respectively. 

Substituting Eq. (14) into Eq. (7) and further into Eq. 
(12) leads to a determined form of the module R n ix,t). 
Using this form we numerically plot the probability den- 
sity R„ versus x for n = 8 at the times t = 0, 7r/2, 2ir 
respectively as Fig. 2a, 2b, and 2c. These plots show 
that the solitonlike train propagates from x c — —10 to 
x c = 10 and approximately keeps its shape. In the next 
half period, 2n < t < Air, an inverse process will occur. 

Collapse and revival of the wave-packets From 
Eq. (13) we know that when small constants A and £>o/co 
are applied to Eq. (12), the wave-packet trains will ap- 
proximately fix their centers. Further we let constant B 
be much greater than A, Eqs. (6) and (7) could lead to 
large change of the function pit) for a — = tt/2. This 
will cause consequently large changes of the widthes and 
highnesses of the wave-packet trains described by Eq. 
(12). The large deformations can be identified with the 
so-called collapse and revival of the wave-packets [23| . 
As an example, we take the parameter set a = 0, [3 = 
-tt/2, A = b = 0.02, B = 10, c = 1, from Eqs. (6), 
(7) and (13) to solve for p{t) and x c (t). The center of the 
wave-packet train could oscillate with amplitude only in 
Abo = 10~ 4 order. The function of width and highness 
pit) periodically changes with amplitude 10 and mini- 
mum value 0.02. The large amplitude and small mini- 
mum value of pit) mean periodical collapse and revival 
of the wave-packet train. Using these parameters, from 
Eq. (12) we plot spatial evolution of the probability den- 
sity with n = 4 for several different times as in Fig 3. 
This figure displays that the initially higher and narrower 
wave-packet train is collapsed to very low and width at 
t = tt and revived to initial state at t = 2ir, as in Fig. 3b 
and 3c respectively. 

In conclusion, we have investigated the quantum and 
classical motions of a time-dependent harmonic oscilla- 
tor, which is associated with a Paul-trapped ion sys- 
tem. A new kind of exact solutions of the quantum- 
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FIG. 2: The probability density of solitonlike train i?| versus x from Eq. (12) for (a) t = 0, (b) t = 7r/2 and (c) t = 2-k. Nine 
wave-packets are shown in propagation with small change of shape. The space-time variables x and t are normalized in units 
lh and cj" 1 , and the probability density is normalized by l^ 1 . 
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FIG. 3: The collapse and revival of wave-packet train from Eq. (12) with n = 4 for (a) t = 0, (b) t = 7T, (c) t = 2k. The 
wave-packet train consisting five packets is shown with approximately fixed center and great deformation. The space-time 
coordinates and probability density are normalized in the same units with Fig. 2. 



mechanical Schrddinger equation is constructed by us- 
ing the corresponding exact solutions of the classical- 
mechanical Mathieu equation, which describes propaga- 
tions and deformations of the wave-packet trains. When 
the deformations are small, they behave like some soli- 
ton trains. The large deformations are identified with 
collapse and revival of the wave-packet trains. If we se- 
lect the initial conditions to fix center of the wave-packet 
trains, the result is reduced to Brown's quantum motion 
Our ground state solution with n = seems to be 



Cessa's coherent state [19(. The state ipi of two wave- 
packets is similar to the Schrodinger's cat state We 
desire the new exact solutions to play an important role 
in treating various harmonically confined systems, e.g. 
the Bose-Einstein condensate held in a magnetic well. 
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